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The /i!/SSM is a supersymmetric model that has been proposed to solve the problems generated 
by other supersymmetric extensions of the standard model of particle physics. Given that R-parity 
is broken in the /ii^SSM, the gravitino is a natural candidate for decaying dark matter since its 
lifetime becomes much longer than the age of the Universe. In this model, gravitino dark matter 
could be detectable through the emission of a monochromatic gamma ray in a two-body decay. We 
study the prospects of the Fermi-LAT telescope to detect such monochromatic lines in 5 years of 
observations of the most massive nearby extragalactic objects. The dark matter halo around the 
Virgo galaxy cluster is selected as a reference case, since it is associated to a particularly high signal- 
to-noise ratio and is located in a region scarcely affected by the astrophysical diffuse emission from 
the galactic plane. The simulation of both signal and background gamma-ray events is carried out 
with the Fermi Science Tools, and the dark matter distribution around Virgo is taken from a A''-body 
simulation of the nearby extragalactic Universe, with constrained initial conditions provided by the 
CLUES project. We find that a gravitino with a mass range of 0.6-2 GeV, and with a lifetime range 
of about 3 X 10^^-2 x 10^* s would be detectable by the Fermi-LAT with a signal-to-noise ratio 
larger than 3. We also obtain that gravitino masses larger than about 4 GeV are already excluded 
in the /xi^SSM by Fermi-LAT data of the galactic halo. 

I. INTRODUCTION 

Evidences indicating the presence of dark matter (DM) can be obtained at very different scales: from cosmological 
ones through the analysis of the angular anisotropies in the cosmic microwave background radiation down to 
galactic scales considering lensing and galaxy dynamics studies. However, all these evidences are not able to provide 
us with complete information on the nature of DM, beyond the fact that it has to be mainly non-baryonic. Since, 
within the standard model of particle physics there are no viable non-baryonic candidates, the existence of DM 
represents one of the most compelling evidences for physics beyond the standard model [2] . 

The "/i from z^" Supersymmetric Standard Model (/ii^SSM) was proposed in the literature to solve the so called 
/i-problem and explain the origin of neutrino masses by simply introducing right-handed neutrinos [3l |4] . It is an 
interesting model that can be tested at the LHC and, as a consequence, its phenomenology has been analized in detail 
recently [nHSj. In the /x;/SSM, i?-parity is broken and therefore the lightest supersymmetric particle (LSP) decays. 
Thus, neutralinos [2] or sneutrinos [7], having very short lifetimes, are no longer viable candidates for the DM of the 
Universe. Nevertheless, if the role of the LSP is played by the gravitino, ^'3/2, its decay is suppressed both by the 
feebleness of the gravitational interaction and by the small i?-parity violating coupling. As a consequence, its lifetime 
can be much longer than the age of Universe and the /iz/SSM gravitino can represent a good DM candidate (51. 

Since the gravitino decays producing a monochromatic photon with an energy equal to half of the gravitino mass, 
its presence can, in principle, be inferred indirectly from the data of the Fermi gamma- ray space telescope [6j . Fermi 
was launched on June 11th, 2008, and its main instrument, the Large Area Telescope (LAT), covers an energy range 
from roughly 20 MeV to 300 GeV with an angular resolution of ~ 0.15 degrees at 68% containment above 10 GeV 
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In the DM halo of the Milky Way the gamma-ray flux coming from DM decay is maximized in the direction of the 
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galactic center where the DM density is larger. This region, however, should be considered with particular care since 
the gamma-ray emission due to the presence of conventional astrophysical sources is not fully understood 

The expected diffuse gamma-ray emission from DM decay in the mid-latitude range (10° < |6| < 20°) was computed 
for a Navarro-Frcnk- White (NFW) profile [TU] in Ref. [5] and compared with the 5-month measurement reported by 
Fermi-LAT The non-observation of sharp monochromatic lines in the gamma-ray spectrum permitted to draw 
bounds on the parameter space of the ^j/SSM gravitino. In particular, values of gravitino mass 7713/2 larger than 
about 10 GeV were excluded, as well as lifetimes T3/2 smaller than about 3 to 5x10^*" s. Notice that because of this 
upper bound on three body decay modes of the gravitino (T^l [T3] are not relevant, and therefore we do not 

considered them in this work. 

It is worth noticing here that in Ref. [14] ([15]) the Fermi-LAT collaboration presented constraints on monochromatic 
emission using 11 (23) months of data for |&| > 10° plus a 20° x 20° square around the galactic center. However, 
the derived limits only refer to the emission above 30 (7) GeV, covering, in the context of the /iJ^SSM, gravitinos 
with masses larger than 60 (14) GeV, leaving our region of interest unconstrained. On the other hand, in work [16) . 
two-years Fermi-LAT data for |6| > 10° have been used to constrain the DM gamma-ray line flux in the energy range 
between 1 and 300 GeV. Lower bounds on T3/2 of about 5 x 10^^ s were obtained in our region of interest below 10 
GeV. Recently, these bounds together with those obtained in [TT] by analyzing the data from EGRET were used in 
[l8] to constrain the parameter space of gravitino dark matter in the bilinear i?-parity violating model. 

We know from X-ray and lensing studies [THl HO] that the main mass component of structures like galaxy clusters 
is DM. Such objects are, therefore, promising targets for DM searches. The main goal of this paper is to explore 
the prospects for detecting /ii^SSM gravitino DM in galaxy clusters using the prediction for 5 years of operation of 
the Fermi-LAT telescope, and taking into account the above mentioned bounds. Although in a very recent work |21j 
based on three years of Fermi-LAT gamma-ray data, the flux coming from nearby clusters was analyzed, the relevant 
bounds obtained on the lifetime correspond to gravitino masses larger than 250 GeV. 

In our analysis we will use the following strategy. The DM density fleld of the nearby extragalactic Universe is 
described using the maps provided in Ref. [22], and based on a constrained iV-body simulation provided by the 
CLUES project [53]. This density is then taken as an input for the Fermi observation simulation tool to predict the 
photon signal. We will use the most recent version of the public Fermi Science Tools [24 to describe the performance 
of the telescope and to simulate both the DM signal and the astrophysical background. From the analysis of the 
simulated photon maps we finally compute our prospects of detection for the /iz^SSM gravitino. 

The paper is organized as follows. In Section [llj the flux of gamma rays from /it/SSM gravitino decay is discussed. 
In Section [Tin the cosmological simulation used to infer the DM distribution in the Local Universe is briefly described, 
while Section |IV| is devoted to summarize the basic of the simulation technique implemented in the Fermi Science 
Tools. Finally, the prospects for DM detection are derived in Section |V| and the conclusions are left for Section \VJ\ 



II. GAMMA-RAYS FROM GRAVITINO DECAY IN THE ^luSSM 

In the supergravity Lagrangian an interaction term is predicted between the gravitino, the field strength for the 
photon, and the photino. Since, due to the breaking of R-parity, the photino and the left-handed neutrinos are mixed, 
the gravitino will be able to decay through the interaction term into a photon and a neutrino [25 . The gravitino 
lifetime T3/2 results to be: 



where |Ky,/|^ is the photino content of the neutrino, and is constrained to be [C/^i/l^ ~ 10""'^^ — 10~^^ in the /xi^SSM, in 
order to reproduce neutrino masses [6 . As a consequence, the gravitino will be very long lived. Additionally, adjusting 
the reheating temperature one can reproduce the correct relic density for each possible value of the gravitino mass 
(see [6] and references therein). 

The detection of DM in several R-parity breaking scenarios has been studied in the literature [S] [551 [55] considering 
the case of gravitinos emitting gamma-rays when decaying in i) the smooth galactic halo, and ii) extragalactic regions 
at cosmological distances. As mentioned in the Introduction, in this work we will analyze the detection of DM in the 
/ij/SSM considering the contribution of gamma-rays coming from in) nearby extragalactic structures. 

In i), the gamma-ray signal is an anisotropic sharp line and the flux is given by 

dE 47rT3/2TO3/2 
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where the halo DM density is integrated along the line of sight /, and we will use a NFW density profile for the Milky 
Way halo compatible with the latest observational constraints as modeled in [27] . Let us remark, nevertheless, that 
since the region that we will study below does not include the inner galaxy, any density profile will give rise to similar 
results. 

On the other hand, in ii), the photons produced by gravitinos decaying at cosmological distances are red-shifted 
during their journey to the observer, and we obtain the isotropic extragalactic flux applying the analysis of Refs. [251128] 
to the iM/SSM. As can be seen e.g. in Figs. 3 and 4 of Ref. p], the sharp line produced by the galactic halo dominates 
over this extragalactic signal. 

Finally, for nearby extragalactic structures in), the gamma-ray signal is a monochromatic line similarly to i), and 
Eq. ([2]) can also be used for the computation of the flux. Actually the contribution from the smooth galactic DM halo 
is practically isotropic in the region around a particular nearby extragalactic structure (at least at high latitudes) 
and less important than the contribution from the object itself. Moreover at high latitudes, the galactic foreground 
is smaller than near the galactic center, so that objects with a lower DM-induced gamma-ray flux can potentially be 
associated to larger prospects of detection than the region near the galactic center. Thus, the study of the extragalactic 
density field is something worth carrying out. In our work this is taken from the results of Ref. [22', where an A^-body 
simulation with constrained initial condition was used. We will review the most important features of this simulation 
in the next Section. 



III. CLUES SIMULATIONS 

CLUES (Constrained Local UniversE Simulations) [23j A^-body simulations aim at describing the formation and 
evolution of DM halos in a way to reproduce, as precise as possible, our Local Universe. To this goal, constrained 
initial conditions are set up using the information from radial and peculiar velocities of galaxies from astrophysical 
catalogs, together with the determination of the masses of the galaxy clusters detected in X-rays [2^ 130] . 

In Ref. [22], a CLUES simulation was used to compute projected Cartesian grid all-sky maps of the DM distribution 
in the Local Universe. In particular, the Boxl60CR simulation was used, containing 1024'^ particles in a box with a 
side of 160 [111[5U]. The characteristics of the most massive clusters such as Virgo, Coma and Perseus, together 
with the Great Attractor, are well reproduced compared to the real objects, apart from a typical mismatch around 
5 Mpc in their position. The all-sky maps are available at Ref. [3T] both in the case of an annihilating and a 
decaying DM particle. Here we use the map corresponding to the case of decaying DM. The gamma-ray flux can be 
derived from the values in the map simply multiplying them by the particle physics factor shown outside the integral 
in Eq. ([2]). This result is used as input for the Fermi- LAT observation simulations, which we will describe in the next 
section. 



IV. SIMULATIONS WITH THE FERMI SCIENCE TOOLS 

For this work, the simulation of gamma-ray events was carried out with the gtobssim routine, part of the Fermi 
Science Tools package v9r23pl. Its output is a list of mock gamma-ray events with corresponding spatial direction, 
arrival time and energy, distributed according to an input source model. Our source model accounts for the gamma-ray 
signal from /ij/SSM gravitino decay as described in previous sections, and for the galactic and extragalactic background 
diffuse emission. In particular, the galactic background emission is mainly correlated with structures in the Milky 
Way since it arises from the interaction of high-energy cosmic rays with the interstellar medium and the interstellar 
radiation field. The far extragalactic background, on the contrary, is supposed to be almost isotropic. Its value is 
based on the modelization of the galactic component, on detected Fermi-LAT sources, and on the solar gamma-ray 
emission. We used the so-called RING model 32 as recommended by the Fermi-LAT collaboration, which is obtained 
as a fit to the real Fermi-LAT data. 

The simulation of gamma-ray events through gtobssim is based on in-flight Instrument Response Functions (IRFs), 
accounting for the telescope effective area, energy dispersion and point-spread function (PSF). Two IRFs publicly 
available are called P6_V3_DIFFUSE and P6_V3_DATACLEAN. They both profit from the improvement in the knowl- 
edge of the telescope performances after the first two years of data taking f33j . The main difference between the two 
IRFs is the fact that DATACLEAN event selection perform most stringent cuts than DIFFUSE on the interpretation 
of an event as a real photon. As will be discussed below, in this work we present results from P6_V3_DIFFUSE, but 
the sensitivity of our results on the choice of the IRFs have been checked without finding any significant effect. 
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FIG. 1; 5'/A'^ all-sky map of gamma-ray emission from gravitino DM decay simulated using the Fermi Science Tools. DM 
events are from the decay of a /i!/SSM gravitino with 7713/2 = 8 GeV and ~ 2.5 x 10^'' s. The signal comes from all three 
DM components: galactic halo, and both, local (from CLUES simulation) and far (isotropic) extragalactic DM density. The 
black square indicates the position of the Virgo cluster. Pixels in the map have an angular dimension of 3 x 3 degrees. 



V. Ai;/SSM GRAVITINO DARK MATTER: PROSPECTS OF DETECTION WITH FERMI-LAT 



Let us finally study the prospects for /it/SSM gravitino DM detection, taking into account the contributions discussed 
in Sect.|ll) First, we simulate the all-sky map of 5-years Fermi-LAT observations, separately for DM gamma-ray events 
and background events. For DM, as an example we show in Fig. |T]the case of a gravitino with mass 7773/2 = 8 GeV 
and lifetime = 2.5 x 10^^ s. Gamma-rays in the energy range between 3.4 GeV and 4.6 GeV are simulated. This 
energy range corresponds to an interval of ±2 AE around the position of the line (4 GeV), where A£'=0.3 GeV is 
the energy resolution (at 68% containment normal incidence) of the telescope at that energy. 

Since we want to determine the Fcrmi-LAT capability of detecting local extragalactic DM structures, the gamma- 
ray emission from DM is what we will refer to as signal. For each direction of the sky tjj, the number of signal (S) 
and background {B) photons are determined integrating over a 3 x 3 degrees region centered in tp. The S/N ratio 
is then defined as S/y/S + B. From the map we can infer which extragalactic structure is the best target to derive 
our prospects for detecting gravitino DM. Let us remark first that we are not simulating the contribution of the point 
sources already detected by Fermi-LAT and present in the 2 year catalog [Mj. Since this contribution represents a 
source of uncertainty, the S/N ratio should not be considered as a good estimator of the prospects for detecting DM 
in the regions where the contamination of the point sources is dominant. Also the discover in the Fermi-LAT data 
of lobes structures |35| , extending up to 50 degrees above and below the galactic plane introduce a new source of 
background not included in this work. The inner galaxy region is characterized by a quite large S/N ratio (around 3). 
However, contamination from point sources as well as the poorly-understood Fermi lobes is known to be large at low 
galactic latitudes, significantly affecting the reliability of our predictions in this region. It is worth noticing here that, 
as mentioned in the Introduction, in [T?! [TS] the region |6| > 10° plus a 20° x 20° square around the galactic center 
was used as target. Let us remark that these works follow a different kind of analysis, searching spectral deviations 
from a power-law behavior. We note that in such a large region of the sky the energy spectra used as background can 
be well described as a power law. Therefore, these works have no problems concerning background uncertainties. 

Thus we neglect the zone with large S/N close to the galactic center and plane, and focus only on the case of 
extragalactic sources. Among these, the object with the largest S/N ratio {S/N — 2.1) is the Virgo cluster. In our 
analysis we only select this cluster, indicated with a black square in Figjl] Although in filamentary regions of the 
cosmic web were pointed out as good targets for DM decay searches, in this analysis, where we are using a different 
particle physics model (and energy range), we do not find in principle a significant S/N ratio in those regions. 

As described in Ref. [22], the projected map is generated with Virgo being fixed in its real observed position because 
it is the best constrained object in the simulation, and therefore it is certainly the best object to consider for our 
purposes. Besides, there are two point sources detected near Virgo, M87 and 2E1228+1437, with an integrated flux 
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FIG. 2: Constraints on lifetime versus mass for gravitino DM in the /i!/SSM. Blue (green) points indicate values of and 
7713/2 of the /ii^SSM gravitino corresponding to a detection of gamma-rays with a. S/N — 5 (3) in the 5x5 degree region centered 
on the position of the Virgo cluster, for a 5 years simulation using the Fermi Science Tools. The blue (green) region indicates 
points with S/N larger than 5 (3). The red dot-dashed line indicates the lower limit on T3/2 obtained from the Fermi-LAT 
measurements of the mid-latitude gamma-ray diffuse emission after 5 months [B]. The yellow dashed line indicates the lower 
limit on r3/2 obtained from the EGRET measurements of the galactic center gamma-ray emission. The black dots show the 
lower limit on r3/2 obtained in the adopted energy bands [IF', from the Fermi-LAT measurements of the |6| > 10° gamma-ray 
diffuse emission after 2 years. The black dashed lines correspond to the predictions of the /ii^SSM 6' for several representative 
values |^7.y^p = 10"^^ 10"^^ 5x 10"", 10"^^ (see Eq. The magenta shaded region is excluded by gamma-ray observations 
such as SPI, COMPTEL and EGRET 1^. 

of (3.3 ± 0.8) X 10^1" ph/cm^s and (1.6 ± 0.6) x IQ-^^ ph/cm^s between 3 and 10 GeV, respectively [M]. In a 5 x 5 
degree region around Virgo and free of the two point sources, the total gamma-ray flux from DM decay is 2.5 x 10"^ 
ph/cm^s (for an example of a 8 GeV gravitino with a lifetime of 5 x lO^^s), one order of magnitude larger than the 
contribution of the point sources. Smaller regions like 1 x 1 or 3 x 3 degrees are more affected by the emission due 
to the point sources. On the other hand, going to larger regions like 7x7 degrees, the signal contribution is not 
significantly increasing and therefore the S/N ratio decreases. Thus, hereafter we select a 5 x 5 degree region around 
Virgo as the best target to obtain our predictions. 

With the purpose of scanning the most interesting portion of the /ii/SSM parameter space, we re-simulate the 
gamma-ray events from the region of 5x5 degrees around Virgo changing the value of the gravitino mass. We run 
17 different simulations of this region, each one with a different value for the gravitino mass, ranging from 0.6 to 10 
GeV, for a given decay lifetime. The lower bound of 0.6 GeV on the gravitino mass is chosen because it corresponds 
to a line energy of about 0.3 GeV, where the PSF (point spread function) of the Fermi LAT becomes larger than 
our region of interest. The energy interval covered by each simulation is [(7713/2/2 — A_E), (7773/2/2 + AE)], where 
the energy resolution (at 68% containment normal incidence) AE is computed at the position of the line. Using the 
results of those simulations, we determine the values of lifetimes corresponding to a S/N = 5 {S/N = 3) . These are 
plotted as blue (green) dots in Fig. [2j as a function of 7773/2. The blue (green) region indicates points with S/N > 5 
(3) . This is the main result of our work. 

Let us remark that the errors in the figure are obtained propagating the statistical errors on the number of signal 
and background events (assuming Poissonian statistics). Regarding possible systematic errors, we note here that 
the points in Fig. [2] are obtained with the P6_V3_DIFFUSE IRF but the simulations were repeated using both, 
P6_V3_DATACLEAN and P7CLEAN_V6, and the results are found to be compatible. The use of P6_V3_DIFFUSE 
allows us to estimate the systematic error on our limits to be between 5 and 20%. 

As mentioned in the Introduction, in Ref. 6J the area below the red dot-dashed line was disfavored by Fermi-LAT 
data of the diffuse gamma-ray galactic emission in the mid-latitude range 10° < |6| < 20°. In addition, in Ref. |16j . 
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the area below the black dots was also disfavored. From a likelihood analysis focused on the region |&| > 10°, lower 
bounds on r3/2 of about 5 x 10^* s were obtained in our region of interest below 10 GeV. On the other hand, the area 
below the yellow dashed line is disfavored by the bounds obtained in [T7] by analyzing the data from EGRET in the 
galactic center region. In particular, we used the upper limits on the gamma-ray line fluxes obtained in that work 
to constrain the /ii^SSM gravitino lifetime. Finally, points in the magenta shaded region are excluded by gamma-ray 
observations from the galactic center obtained with the SPI spectrometer on INTEGRAL satellite, and the isotropic 
diffuse photon background as determined from SPI, COMPTEL and EGRET data [56] . 

On the other hand, the black dashed lines correspond to the predictions of the /iJ^SSM for several representative 
values of the R-parity mixing parameter. As mentioned in Sect. Hi this is constrained to be jJ/^iyp ^ 10~^^ — 10^^^ 
in the yiti'SSM |6j, in order to reproduce the correct neutrino masses. As a consequence, any acceptable point must 
be in the area between the left and right black dashed lines. Let us remark, however, that these bounds are very 
conservative, as discussed in [6 , and in fact the results of a scan of the low-energy parameter space of the /ii^SSM 
implied that the range 10^^^ < iKyj/P < 5 x 10~^^ is specially favored. The corresponding lines are also shown in the 
figure. 

The combination of the constraints associated to red dot-dashed and black dashed lines, implies already that values 
of the gravitino mass larger than about 10 GeV are excluded, as well as lifetimes smaller than about 3 to 5x10^^ s 
[6]. Actually, in the region of gravitino masses between 0.6 and about 1.5 GeV, lifetimes smaller than about 7 to 3 
X 10^"^ s, respectively, are excluded because of the constraints associated to the yellow dashed line. When constraints 
associated to black dots are also imposed, it turns out that the gravitino mass has to be smaller than about 4 GeV, 
and lifetimes have to be larger than about 6 x 10^^ s for gravitino masses between 2 and 4 GeV. Thus, the combination 
of these results with the one obtained above for detection of DM from Virgo in 5 years of Fermi-LAT observations, 
leaves us with the blue and green areas above the yellow dashed and red dot-dashed lines, and gravitino mass smaller 
than 2 GeV, as those with good prospects for DM detection. Summarizing, we find that a gravitino DM with a mass 
range of 0.6-2 GeV, and with a lifetime range of about 3 x 10^^-2 x 10^^ s would be detectable by the Fermi-LAT with 
a signal-to-noise ratio larger than 3. If no gamma-ray lines are detected in 5 years, these regions of the parameter 
space of the /ii^SSM would be excluded. 

VI. CONCLUSIONS AND OUTLOOK 

In this work we have obtained the regions of the parameter space (m3/2, T3/2) of the /zr^SSM with the best prospects 
for the detection of a gamma-ray monochromatic line from the decay of gravitino DM (see Fig. [2]). Summarizing, 
we find that a gravitino DM with a mass range of 0.6-2 GeV, and with a lifetime range of about 3 x 10^^-2 x 10^* 
s would be detectable by the Fermi-LAT with a signal-to- noise ratio larger than 3. We also obtain that gravitino 
masses larger than about 4 GeV are now disfavored in the /ii/SSM by Fermi-LAT data of the galactic halo. 

In the analysis we have assumed 5 years of observation of the Virgo cluster by the Fermi-LAT space telescope. 
This cluster was selected as our optimal target due to its particularly high S/N ratio. Of course, a more precise 
determination of the Fermi-LAT possibilities of detecting gamma-ray lines towards Virgo would require the simulation 
of both, M87 and 2E1228-I-1437. Also a more sophisticated analysis pipeline than the computation of the S/N ratio, 
possibly through the determination of a test-statistics (TS) likelihood in order to derive the lower value of T3/2 for 
which the signature of a line would be detectable with respect to the background. 

Let us remark that the simulation of the gamma-ray flux was carried out with the use of the gtobssim routine from 
the Fermi Science Tools, whereas the DM distribution around the cluster has been modeled following the results of 
Ref. [22j based on a constrained A^-body simulation from the CLUES project [23j. With the present work we have also 
confirmed the potential of using extragalactic massive structures as optimal targets for decaying DM detection. For 
such a goal the maps of the local extragalactic DM distribution produced in Ref. [22] represent a unique, ready-to-use 
tool. 

We conclude that there are good prospects for Fermi to detect monochromatic lines from gravitino decay in the 
energy range spanning from hundreds MeV to few GeV. That is also the energy range where it is more difficult 
to extract information from the data, due to imperfect parametrization of the background as a simple power law. 
Nevertheless, our results in Fig. [2] can be considered as an additional motivation to extend the Fermi-LAT analysis 
on lines to energies below 2 GeV. 
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